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Using a sample of 152 million BB pairs accumulated with the Belle detector at the KEKB 
e"'"e~ collider, we study the decay mechanism of three-body charmed decay B~ —> A+pvr". The 
intermediate two-body decay B~ Ec(2455)''p is observed for the first time with a branching 
fraction of (3.7 ± 0.7 ± 0.4 ± 1.0) x 10~^ and a statistical significance of 8.4a. We also observe a 
low-mass enhancement in the (A^p) system, which can be parameterized as a Breit-Wigner function 
with a mass of (3.35+°;o2 ± 0.02) GeV/c^ and a width of (0.07to;o3 ± 0.04) GeV/c^ We measure 
its branching fraction to be (3.9lg'7 ± 0.4 ± 1.0) x 10"^ with a statistical significance of 6.2a. The 
errors are statistical, systematic, and that of the pK^Tv^ decay branching fraction. 

PACS numbers: 13.25.Hw, 14.20.Lq 



Recently three-body baryon production in charmless 
B decays has been studied with the BeUe detector 
P) S S Si- Analysis of these decays shows a common 
feature: the invariant mass of the baryon-antibaryon sys- 
tem is peaked near threshold. This feature has generated 
much theoretical discussion and may be due to a frag- 
mentation effect, production of resonances near thresh- 
old or final state interaction of the produced baryon- 
antibaryon system 5, 6, 7, 8, 9). It is of interest to learn 
whether similar behavior is also observed in B decays to 
charmed baryons. The three-body decay B~ A+p7r~ 
has been previously studied at CLEO and Belle [llj 
with 9.2 fb~^ and 29.1 fb^^ of data, respectively. Here we 
report analysis of the B^ A+p7r~ decay on a Dalitz 
plane based on a data sample of 140 fb~^ accumulated 
at the T(4S') resonance with the Belle detector at the 
KEKB e+e- collider ^\. 



The large-solid-angle magnetic spectrometer Belle de- 
scribed in detail elsewhere [l^l consists of a three-layer 
silicon vertex detector (SVD), a 50-layer cylindrical drift 
chamber (CDC), a mosaic of aerogel threshold Cerenkov 
counters (ACC), a barrel- like array of time-of-flight scin- 
tillation counters (TOF), and an array of CsI(Tl) crystals 
(ECL) located inside a superconducting coil providing a 
1.5 T magnetic field. An iron flux return located out- 
side the coil is instrumented to detect muons and 
mesons (KLM). We use a GEANT based Monte Carlo 
(MC) simulation to model the response of the detector 
and determine its acceptance 14 1. 



the combined response of the CDC, ACC and TOF. We 
require the impact parameters of all primary tracks with 
respect to the interaction point (IP) to be within ±1 cm 
in the transverse {x — y) plane and within ±5 cm in the 
z-direction (opposite to the e"*" beam). Proton, kaon and 
pion candidates are selected using p/K/i: likelihood func- 
tions provided by the PID system. We require the like- 
lihood ratios Li/{Li + Lj) to be greater than 0.6, where 
the subscript i denotes the selected particle and j the 
other two particle species. The PID eflaciency is 98% for 
each track, and the fake probability of a pion (kaon) to 
be identified as a kaon (pion) is less than 5%. The proba- 
bility for a pion or kaon to be identified as a proton is less 
than 2%. We detect the A+ via five decay modes: A+ — > 
pK~T:^, pK^, Att+j pif 7r+7r~ and Att^tt+tt". Inclu- 
sion of charge conjugate states is implicit unless otherwise 
stated. Neutral kaons and A baryons are reconstructed 
in the ifg— >7r"'"7r~ and A— >p7r~ decay, respectively. The 
B^ — !■ AjtpTT" events are identified by their energy differ- 
ence AE = i^Ei) — Sbcam, and the beam-energy con- 
strained mass Mbc = V ^Lam ~ (EpO^^ where Ehoam is 
the beam energy, and pi and Ei are the three-momenta 
and energies of the B meson decay products, all defined 
in the e'^e~ center-of-mass system. To suppress contin- 
uum background, we impose requirements on the angle 
between the thrust axis of the B candidate tracks and 
that of the other tracks and on the ratio of the second to 



the zeroth Fox- Wolfram moments [15 1 



The event selection is based on track information from 
the SVD and CDC and particle identification (PID) from 



Fig. [T] shows the AE distribution for Afbc > 
5.27 GeV/c^ for the selected B candidates. The signal 
yield is extracted by a fit to the AE distribution, as it 
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FIG. ] 



AE distribution for A/bc > 5.27GeV/c^ for 
candidates. The curve shows the result of the fit. 



is free from combinatorial backgrounds from other B de- 
cays. The fit uses a Gaussian for the signal fixed to MC 
data and a third order polynomial for the background for 
which a MC study shows a broad background below the 
signal due to continuum events and combinatorial back- 
grounds from other B decays, such as — > A+p7r°. We 
obtain 264 ± 20 signal events with a statistical signifi- 
cance of 17.1 a. The uncertainty due to the background 
parameterization is estimated by repeating the fit with 
first and third order polynomials and found to be small 
(1.5%). The significance is defined as ^y—2\n{Co/jCniax), 
where Cmax and Cq denote the maximum likelihoods with 
the fitted signal yield and the yield fixed at zero, respec- 
tively. 

Fig. [Ha) shows the M(A+7r^) distribution, where 
clear signals are seen from the intermediate two-body 
B- Sc(2455/2520)°p decay [13, [Hi. The open his- 
togram is the distribution from the B signal region 
{\AE\ < 0.03 GeV and Mbc > 5.27GeV/c2). The 
hatched histogram is the distribution from sideband re- 
gions (-0.10 GeV < A£: < -0.04 GeV or 0.04 GeV < 
AE < 0.20 GeV) normalized to the B signal region area. 
The curve shows the result of the fit, which includes 
the contribution from Ec(2455/2520)" A+tt" decays 
and the background parameterized with a linear func- 
tion. The I]c(2455/2520)'^ signal shapes are fixed from 
MC assuming a Breit-Wigner function convolved with 
the resolution function and using Ec(2455/2520)'^ masses 
and widths from Ref. [l3|. To extract the yields for 
i3~—>Ec (2455/2520)°p decays, background from contin- 
uum and/or other B decays is taken into account by fit- 
ting simultaneously the B signal and sideband regions 
in the AE distribution. From the fit, we obtain 35.3lgo 
signal events with a statistical significance of 8.2 a for the 
B~ i;c(2455)°p decay, and 12.6t4;7 signal events with 
a statistical significance of 3.0 a for the B~ Ec(2520)''p 
decay. 

Fig.HKb) shows the M{pn~) distribution for the B^ 
A+p7r~ candidate events from fits to the AE distribu- 
tion in 100 MeV bins of (p7r~) mass with constraints 
of M(A+7r-) > 2.6GeV/c2 to remove E° intermedi- 
ate states and Af(A+p) > 3.5GcV/c^ to remove an en- 
hancement at low (A+p), which is discussed below. The 
histogram shows the result of a fit including the fol- 
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FIG. 2: (a) Af(A+7r") distribution for the B signal re- 
gion (open histogram) and fit results (curve); the distri- 
bution in the sidebands is also shown (hatched). (b,c) 
B~ yields (points) from fits to AE distributions, (b) in 
bins of M{pTT~), requiring MiA+n-) > 2.6 GeV/c^ and 
M(A+p) > 3.5 GeV/c^ and (c) in bins of A/(A+p), requiring 
A^(A+7r~) > 2.6 GeV/c^ and M{p-K-) > 1.6 GcV/c^ The 
histograms show fit results, see the text. 



lowing contributions: three-body phase space, B 
A+A(1232)— [13, and two other contributions, with 
parameters close to the A(1600) and A(2420) resonances 
tentatively referred to as Ax(1600) and Ax(2420). The 
signal shapes are fixed from the MC data. Both three- 
body phase space and A(1232) contributions have a 0.5f7 
significance. The statistical significance of the Ax(1600) 
contribution is 6.7a with a yield of 82 ± 12 events while 
that of the Ax(2420) is 4.7o- with a yield of 41 ± 9 events. 

Fig- [2Ic) shows the M (A+p) distribution for B~ — > 
A+]57r~ decay candidate events from fits to the AE distri- 
bution in 50 MeV bins of (A+p) mass with Af(A+7r~) > 
2.6GeV/c'^ to remove Ec contributions and M{pTr^) > 
1.6 GeV/ to remove low {pTT~) masses. A low 
mass enhancement is observed. The histogram is 
the result of a fit parameterizing the distribution 
with a Breit-Wigner peak and feed-downs from the 
B- A+Ax(1600/2420)"~ contributions. This fit 
gives a mass of {3.35to°ol)GeV/c^ and fuU width of 
(0-07lo;o3) GeV/c2. The fit yield is 50 ± 10 events with 
a statistical significance of 5.6f7. A second peak near 
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3.8 GeV/c^ has a mass of (3.84 ±0.01) GeV/c^ and width 
of (0.03 ± 0.03) GeV/c2. The yield of this peak is 15 ± 6 
events {2.8a) and not studied any further. 

Systematic uncertainties are estimated by performing 
fits with different background parameterizations includ- 
ing the contributions of the B- A+Ax(1600/2420) — 
feed-downs with a free number of events or a broad Breit- 
Wigner function with and without the second peak at 
3.8 GeV/c^. The mass variation is less than 0.02 GeV/c^, 
and the width varies by less than 0.04GeV/c^. 

For the (A+p) structure, we studied the distribution of 
the helicity angle, 9(A+p), defined as the angle between 
the A+ momentum and the direction opposite to the B 
meson momentum in the (A+p) rest frame. If the (A+p) 
structure is due to fragmentation, the distribution will 
be asymmetric, while for a resonance it will be symmet- 
ric y and could provide information on the spin of the 
(A+p) state. Fig. [3] shows the efficiency corrected helic- 
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FIG. 3: The efficiency corrected helicity distribution for the 
(A^p) structure. The solid line is the the result of the fit to 
(1 -^acos^e(A+p)). 

ity distribution for this decay for events from the region 
Af(A+p) < 3.6 GeV/c^, where the data points were ob- 
tained from fits to the AE distributions. The data is 
consistent with the uniform distribution expected for a 
J = state (x^/ndf — 0.97). A fit of our data to a general 
formula for the angular distribution (I -f acos^ 0(A+p)) 
for a J = 1 state [H HH gives a = (-0.15 ± 0.54) 
(x^/ndf — 1.12), see the solid curve in Fig. [31 The ob- 



served helicity asymmetry is 



JV4 



0.32±0.14, where 



N++N- 

and N_ are the efficiency corrected numbers of events 
with cos0(A+p) > and < 0, respectively. No definite 
conclusions can be drawn from the helicity studies with 
this statistics. 

We determine the contributions from Ec (2455/2520) 
and the low (A+p) mass structure, taking into account 
cross-talk between different resonant states and the vari- 
ation of detection efficiency on the Dalitz plane. Fig. [4] 
shows the Dalitz plot of M(p7r~)^ vs M(A+7r~)2 for 
the events in the B signal region of \AE\ < 0.03 GeV 
and Mbc > 5.27GeV/c^, where we estimate a back- 
ground contamination of 37% in total. The Dalitz plot 
is subdivided into six regions corresponding to the six 
states discussed above: I) Ec(2455)°p — M(A+7r") < 
2.48GeV/c2; 2) Ec(2520)"p — M(A+7r-) > 2.48GeV/c2 



and M(A+7r-) < 2.6GeV/c2; 3) A+A(1232)— — 
Af(A+7r-) > 2.6GeV/c2 and M{p7r~) < 1.4GeV/c2; 
4) A+Ax(1600)-~ — M(A+7r-) > 2.6GeV/c2, 
M{Ptt-) > 1.4GeV/c2 and M(p7r-) < 2.0GeV/c2; 5) 
A+Ax(2420)~ — M(A+7r-) > 2.6GeV/c2, M{ptt~) > 
2.0GeV/c2 and M{A+p) > 3.6GeY/c^; 6) (A+p) en- 
hancement — M(p7r") > 2.0GcV/c2 and M{A+p) < 
3.6GeV/c2. 

The resonance parameters are taken from [16]. The 
(A^p) structure is represented as a Breit-Wigner func- 
tion with mass 3.35GeV/c2 and full width 0.07GeV/c2, 
visible as a band in region 6 of Fig. |4l The B signal 
yield in the i-th region is given as Xi = ^ij ' ^ j 

where Sij is the probability to reconstruct the j-th in- 
termediate state in the i-th Dalitz plot region (estimated 
from the MC sample of j-th state). We extract the sig- 
nal Yj from a simultaneous fit of the B signal yields Xi 
for the six AE distributions, where the width is fixed 
from MC data. Then, we calculate the branching frac- 
tion Bj = Yj/{NBB ^i^t pK-TT+)) for the j-th 
intermediate state, combining the B signals tagged with 
the five A+ decay modes, is the combined efficiency 
given by ^r^fe X rfc(A+)/r(A+ -> pK-^+) 0. 77^ is 
the efficiency of the B signal with the fc-th A+ decay 
determined from MC. 
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FIG. 4: The M(p7r")^ vs M(A+7r")^ Dafitz plot subdivided 
into six regions. The points are the data in the B signal 
region. The shaded area is the phase space MC data. 

The resulting yields, branching fractions and statisti- 
cal significances for each intermediate state are listed in 
Table IT] For the branching fractions, the first error is sta- 
tistical and the second is systematic, whereas the third 
(26%) comes from the uncertainty in B{A'^^pK~Tr'^). 

Systematic uncertainties in the detection efficiencies 
arise from the track reconstruction efficiency (5 — 7% de- 
pending on the process, assuming a correlated systematic 
error of about 1% per charged track); the FID efficiency 
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(about 7% assuming a correlated systematic error of 2% 
per proton and 1% per pion or kaon); and MC statistics 
(1 — 2%). The other uncertainties are associated with 
r(A+)/r(A+ ^ pK-iT+) (1 - 2%); the number of BB 
events (0.8%); and the parameters of the low mass (A+p) 
enhancement, which can contribute up to 6% to the un- 
certainty of its branching fraction. The total systematic 
error is estimated to be 9 — 11% depending on the inter- 
mediate state. 

We estimate separately a possible effect of interference 
between the different observed intermediate states. Using 
special MC samples in which each of the relative phases 
among the six states is varied in steps of 90°, we com- 
pare the signal yield in individual regions of the Dalitz 
plot to that obtained using simulated events without any 
interference. The maximum deviation is treated as the 
uncertainty due to interference and is given in Table [J 
This simplified treatment of interference does not take 
into account a possibility of reduced compatibility be- 
tween the simulated distributions with interference and 
data and indicates that the low mass (A+p) enhancement 
can be partially described by such an effect. 

In summary, using a sample of 152 million BB pairs, 
accumulated with the Belle detector at the KEKB col- 
lider, we performed an analysis on the Dalitz plane of 
the three-body charmed decay B^ A+p7r~. We 
report first observation of the two-body decay mode 
B~ — > I]c(2455)°p and measure its branching fraction 
to be (3.7 ± 0.7 ± 0.4 ± 1.0) x 10"^ with a statistical sig- 
nificance of 8.4cr. We also observe a low mass enhance- 
ment in the (A^p) system, which can be parameterized 
as a Breit-Wigner function with a mass of (3.35j;Q Q2 ± 
0.02) GeV/c2 and a width of (0.07t^;^3 ± 0.04) GeV/c^. 
The branching fraction of the B~ decay to this structure 
is (3.9lo 7 ± 0.4 ± 1.0) X lO"'^ with a statistical signifi- [15 



cance of 6.2a. The current data are not sufhcient to de- 
termine an origin of this enhancement. The total three- 
body B~ — > A'^p'K~ decay branching fraction has been 
measured to be (20.1 ± 1.5 ± 2.0 ± 5_^x 10-^ which is 
consistent with previous results 1^, 11 1. The branching 
fractions measurements supersede those in Ref. [ll| . 
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TABLE I: Branching fractions for intermediate two-body states in B — >A+p7r decay. 
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